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ABSTRACT 

We compare the WMAP temperature power spectrum and SNIa data to models with 
a generalized Chaplygin gas as dark energy. The generalized Chaplygin gas is a com- 
ponent with an exotic equation of state, px = —A/p x (a polytropic gas with negative 
constant and exponent). Our main result is that, restricting to a flat universe and 
to adiabatic pressure perturbations for the generalized Chaplygin gas, the constraints 
at 95% CL to the present equation of state wx — Px I Px and to the parameter a 
are — 1 < wx < —0.8, < a < 0.2, respectively. Moreover, we show that a Chaply- 
gin gas (a = 1) as a candidate for dark energy is ruled out by our analysis at more 
than the 99.99% CL. A generalized Chaplygin gas as a unified dark matter candidate 
(flcDM = 0) appears much less likely than as a dark energy model, although its x 2 is 
only two sigma away from the expected value. 
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1 INTRODUCTION 

The picture of a universe mostly filled by an unknown component, Dark Energy (DE henceforth), seems overwhelmingly 
justified by the latest observations of CMB anisotropics (Spergel et al. 2003). 

While there is a strong degeneracy of the multitude of DE candidates for the background evolution - and therefore in the 
SNIa data (Perlmutter et al. 1999; Riess et al. 1998) -, additional information should be contained in the linear fluctuations. 
However, in order to differentiate among DE models, the avenue of studying DE fluctuations is not promising as it seems. 
Infact, in the simplest uncoupled quintessence models (Ratra & Peebles 1988; Frieman et al. 1995; Caldwell, Dave & Steinhardt 
1995) based on scalar fields with ordinary Klein-Gordon Lagrangian, quintessence fluctuations are evanescent on sub-horizon 
scales. This is due to the presence of a Jeans scale for scalar field perturbations, of the order of the Hubble radius. This Jeans 
scale is related to a unitary sound speed (in c = 1 units) for scalar field perturbations. For the same reason, isocurvature 
perturbations caused by an offset of the quintessence component (Abramo & Finelli 2001) are not observationally harmful as 
other isocurvature modes. 

A possibility of having a smaller Jeans scale for dark energy fluctuations is offered by theories with non standard kinetic 
terms for the scalar field candidate. A scalar field with a Born-Infeld action is an example: 

s = J d ^vf^^, (1) 

where M is a fundamental scale and V is a potential. Recently, the tachyon field (with V = V(ip)) motivated from string 
theory has received a lot of attention (Sen 2002, Gibbons 2002, Padmanabhan 2002). A Chaplygin gas is a hydrodinamical 
description of the scalar field with the Born-Infeld Lagrangian in Eq. (1) with constant V (Kamenshchik, Moschella & Pasquier 
2001). A Chaplygin gas is characterized by a pressure which is inversely proportional to the energy p oc — 1/p. A generalized 
Chaplygin gas (henceforth GCG) is a perfect fluid with a polytropic equation of state (Bento, Bertolami & Sen 2002): 
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where a is a parameter between and 1 in order to have a sound speed at most luminal for perturbations. A is a constant with 
dimensions [M 4(1+q) ]. A GCG is a phenomenological extension of the Chaplygin gas, which exhausts all the possibilities for 
a polytropic perfect fluid dark energy candidate, whose perturbations are stable on small scales (Fabris, Goncalves & Souza 
2002a, Carturan & Finelli 2002 - henceforth CF -). 

The behaviour in time of a GCG interpolates between dust and a cosmological constant, with an intermediate behaviour 
as p — ap. The Jeans instability of GCG perturbations is first similar to CDM fluctuations (when the GCG has a negligible 
pressure) and then disappears (when the GCG behaves as a cosmological constant) (see CF). Both this late suppression of 
GCG fluctuations and the apperance of a non zero Jeans length leave a large integrated Sachs- Wolfe (ISW) imprint on the 
CMB anisotropics (CF). A model in which the GCG is the only dark component is referred to as a unified dark matter/dark 
energy (UDM) model. 

The viability of GCG as DE and UDM model has been first analyzed in the context of Supernovae (Fabris et al. 2002b; 
Avelino et al. 2002; Makler, de Oliveira & Waga 2003) and other complementary data (Dev, Jain & Alcaniz 2003, Alcaniz, 
Jain and Dev 2003). The comparison with CMB data has been the next step (CF, Bean & Dore' 2003) and it has led to much 
stronger constraints (Bean & Dore' 2003). 

The imprint of a GCG is also present on the matter power spectrum, since GCG perturbations affect both CMB 
anisotropics and structure formation. This last issue has been considered by Sandvik et al. (2002) in the context of UDM 
models, claiming that a GCG is ruled out as a UDM candidate at 99.999% (Sandvik et al. 2002). However, their analysis 
does not take into account baryons, while we will show they are relevant for the shape of the total matter power spectrum. 
Moreover, the analysis by Sandvik et al. is based on a linear treatment of perturbations until the present time, neglecting 
any non linear effects which may be important and unexpected because of the time dependence of the GCG Jeans length. 
Such nonlinear effects should be more important for the present matter power spectrum than for the CMB spectrum. For this 
reason, we focus on the comparison with CMB data. 

In this paper we compare the predictions of GCG models with the recently released WMAP data, which provide the most 
precise determination of the CMB spectrum so far (Hinshaw et al. 2003), and with SNIa data (Perlmutter et al. 1999, Riess 
et al. 1998). 



2 THE MODEL 

In this section we review a brief and basic description of the model (for more details see CF). The conservation equation for 
the GCG in a Robertson- Walker metric is solved by 

^ = ^ + -^1)]^' ( 3 ) 

where a is the scale factor (a = 1 today) and B is an integration constant with dimension [M 4( - 1+a - ) ]. 
The equation of state is therefore 

/ \ Px A 

wx(a) = -=--^, (4) 

while the present equation of state is 
A 

Wx = -aTb- 

Instead of A, B we find convenient to employ the parameters fix , wx , where fix is the present value of the density parameter 
of the GCG and 

a = - Wx (n x pc) 1+a , 

B = (l + w x )(tl x p c ) 1+a , 

where p c is the present critical density. Beside the GCG we introduce for generality also baryons and CDM. Although the 
GCG is very interesting as an unified model of dark energy and dark matter, it is conceivable that it is in fact only a DE 
component additional to ordinary CDM and baryons. 

Perturbations of GCG are stable on small scales since the sound speed is positive (c x = dpx/dpx = —awx). An 
important consequence of assuming a perfect fluid instead of a scalar field as a candidate of dark energy is the absence of 
intrinsic non adiabatic pressure perturbations. Infact, for a GCG pressure perturbations are locked to density perturbations 

Spx = c x $px ■ (5) 
Instead, in presence of non adiabatic pressure perturbations, we have in general: 
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Figure 1. The Ci spectrum varying wx is plotted in comparison with a ACDM model (the upper curve): wx — 
—0.8, —0.91, —0.99, —0.999 from bottom to top, respectively. The other parameters are a = 1, h = 0.72, w b = 0.024, w c = 0.14, 
n s = 1. The spectra are COBE normalized at t = 10. 



Spx = c x 5px + Apx ■ (6) 

For example, in a scalar field model with standard kinetic term for quintessence <? x = 1 and Apx = 2Vtf,5(f>, where V4, is the 
first derivative of the potential and S<j> is the field fluctuation. 

A perfect fluid description of a GCG offers the chance to explore the possibility of the basic relation (5) for dark energy 
pressure perturbations. 

In Figs. 1 and 2 we compare the Ci temperature spectrum with the WMAP data varying wx and a, respectively (see 
also CF). These figures show how strongly the Ci spectrum depends on perturbations of GCG, when this is a DE candidate 
(Q c is fixed to 0.27). Fig. 1 shows how the resulting CMB spectrum strongly differ even for models which have wx < —0.91. 
Notice that the 2nd and 3rd peak move to the left increasing wx, because then the expansion is less and less accelerated, 
while the 1st peak remains more or less stable because of the concurring strong integrated Sachs-Wolfe (ISW) effect. Fig. 2 
shows the strong dependence of CMB spectrum on a, which regulates the Jeans length of GCG perturbations and therefore 
the amount of ISW effect. The peak-to-plateau ratio decreases for a increasing, because after the dust behaviour the GCG 
induces a strong ISW. The dependence of CMB spectrum and LSS on the Jeans length of DE perturbations is a generic effect 
and it has been also found recently in K-essence models (DeDeo, Caldwell & Steinhardt 2003), but it is less severe than in 
GCG models. 



3 LIKELIHOOD ANALYSIS 

We compare the models to the combined power spectrum estimated by WMAP (Hinshaw et al. 2003). To derive the likelihood 
we adopt a version of the routine described in Verde et al. (2003), which takes into account all the relevant experimental 
properties (calibration, beam uncertainties, window functions, etc). Since the likelihood routine employs approximations that 
work only for spectra not too far from the data, we run it only for models whose x 2 P er d.o.f. is less than 1.5 from the number 
of degrees of freedom. 
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Figure 2. The Cg spectrum varying a is plotted in comparison with a ACDM model (the upper solid line): a = 1, 0.5, 0.2, from bottom 
to top, respectively. The other parameters are wx — —0.99, h = 0.72, wi, = 0.024, w c = 0.14, n s = 1. The spectra are COBE normalized 
at I = 10. We note how the ACDM curve and a = are very close, but not identical because of the different perturbation sectors. 



Our theoretical model depends on two GCG parameters, four cosmological parameters and the overall normalization A^: 

a,w x ,n s ,h,uj b ,uj c ,N. (7) 

where u>i,,c = Qb,ch 2 and n s is the primordial fluctuation slope. As anticipated, to save computing time we found necessary to 
restrict the analysis to a flat space; moreover, we fixed the optical depth to t — 0.17, the mean value found by WMAP (Spergel 
et al. 2003). The overall normalization has been integrated out numerically. We calculate the theoretical Ce,t spectra by a 
modified parallelized CMBFAST (Seljak & Zaldarriaga 1996) code that includes the full set of perturbation equations (see CF). 
We do not include gravitational waves and the other parameters are set as follows: T cm b = 2.726K, Yne ~ 0.24, N„ = 3.04. 

We evaluated the likelihood on a non-uniformly spaced grid of roughly 50, 000 models (for each normalization) with the 
following top-hat broad priors: wx <E (-1,-0.5), a 6 (0, .5), n s <E (0.8, 1.2), u>b £ (0.005, 0.04), lu c e (0.0, 0.2). For 
the Hubble constant we adopted the top- hat prior h £ (0.5,0.9); we also employed the HST result (Freedman et al. 2001) 
h = 0.72 ±0.08 (Gaussian prior). One reason to adopt a grid approach than a Markov chain method (as in Verde et al. 2003) 
is that the independent grid evaluation can be parallelized with maximum efficiency. 



4 RESULTS FROM COMPARISON WITH WMAP 

In Fig. 3 we plot the likelihood functions for each parameter, marginalizing in turn over the others. The horizontal lines mark 
the 68% and 95% CL; the vertical lines the 68% and 95% upper bounds. We obtain the bounds 

a < 0.2(0.05) 

w x < -0.8(-0.92) (8) 

at 95%(68%) CL. These constraints are the most stringent ones obtained so far using only CMB data; the limit on a is 5-10 
stronger than the previous one obtained with pre- WMAP data (Bean & Dore' 2002). For the other parameters we obtain 
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Figure 3. Marginalized likelihood functions for the six cosmological parameters. The dotted lines are for the HST prior on the Hubble 
constant. The horizontal long-dashed lines are the confidence levels at 68% and 95%. The vertical long-dashed lines in the panels for 
a , wx mark the upper bounds at 68% and 95% CL. 



results not very different than the standard ones. We also show the results imposing the HST prior h — 0.72 ± 0.08; they are 
almost identical to the one without prior on h. The best fit DE model is n s = 1, a = 0, wx = —0.98, h — 0.71, uib = 0.024, 



= 0.12 has x /d.o.f.=979/892. 



5 CONSTRAINTS ON THE UNIFIED DARK MODEL 

Any DE model which interpolates in time between a dust epoch and an accelerating stage has the possibility of playing the 
role of a UDM candidate, if one renounces to the presence of CDM. Naturally, in such a case one should explain also why 30% 
or so of the matter collapsed into galaxies and clusters while the remaining did not. The GCG has received a lot of attention 
also in this light of a UDM model. While SNIa data allow for this possibility, it was claimed that LSS data completely ruled 
it out (Sandvik et al. 2002). Here we find that WMAP data alone (for which a linear analysis is sufficient) show that a 
GCG as a unified candidate is much less likely than a model which includes a CDM component. The likelihood for lo c = is 
infact roughly 50 times smaller than for ui c = 0.12; formally, this rules out the UDM at slightly more that 99.99% C.L. when 
compared with GCG as DE. The lowest value allowed by WMAP data is 0.1 (95C.L.). However, if we calculate the \ 2 f° r the 
best fit among the UDM cases (corresponding to a model with n s — 0.98, a = 0, wx = —0.8, h = 0.82, cut = 0.026) we find 
X 2 /d.o.f.=986/893, indicating that the UDM model is an acceptable fit to WMAP data at at 2a level. 

Note that our results have been obtained varying also n s , which has been instead fixed in the analysis of Bean & Dore' 
(2002) and of Sandvik et al. (2002) (who also fixed h). 
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Figure 4. Likelihood in the plane wx,ct. The contours are at the 68%, 95% and 99% CL, inside to outside. 



6 ADDING THE SUPERNOVAE CONSTRAINTS 

We now proceed with a further refinement of our result, by taking into account the constraint from the SN I a Hubble diagram. 
This has been evaluated by several groups (Fabris et al. 2002a; Avelino et al. 2002; Makler et al. 2003), but not in a suitable 
form for our purposes. 

The luminosity distance in a GCG model in flat space is 



d L (z) = (l + z) 



r dz' 
J H{ z >y 



where the Hubble function is 



H 2 (z) = Hi 



+ zf + (1 - Q m )((l + wx)(l + z) 3{1+a) - w x )rh 



and Q. m — fit + Q c . 

We compare the luminosity distance to the SCP data of Perlmutter et al. (1999) (their fit C), to which we add the 
supernova at z w 1.7 (Benitez et al. 2002). We show in Fig. 6 the 2-dimensional likelihood function marginalized over the 
"nuisance parameters" (see definitions in Perlmutter et al. 1999) a, M* , and over fi m with Gaussian prior Q m = 0.3 ± 0.1. 
Notice that the contours are almost insensitive to a since (z) at small redshifts is independent of a. This shows how CMB 
data constrain the GCG models much more than SNIa data. In Fig. 7 we multiply the SN and WMAP likelihood functions. 
The final constraints turn out to be almost identical to those in (8). 



7 CONSIDERATIONS ON THE MASS POWER SPECTRUM 

The influence of the GCG on structure formation is another important observational imprint of this candidate for DE and 
UDM model. GCG fluctuations accumulate on all scales when the GCG behaves as CDM. As soon as the GCG equation 
of state drops below zero, the GCG fluctuations are no longer Jeans unstable at all scales, oscillating and decaying. This 
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Figure 5. Likelihood functions for the GCG model (full curve) and the unified model (Q c = ; dashed curve). The horizontal long- 
dashed lines are the confidence levels at 68% and 95%. The vertical long-dashed lines in the panels for a , wx mark the upper bounds at 
68% and 95% CL. 



exotic behaviour of fluctuations occurs also when the GCG simulates A. As already said in the introduction, this suppression 
of fluctuations is in part responsible for the large ISW effect in CMB anisotropies. It is clear that a similar behaviour of 
fluctuations differs both from CDM and standard quintessence, and it is more important in absence of CDM, i.e. in a UDM 
model. 

Sandvik et al. (2002) used a comparison of the predicted power spectrum in GCG models with LSS data - 2dF survey 
(Colless et al. 2001; see also Tegmark, Hamilton & Xu 2002) - to rule out the GCG as a UDM candidate at 99.999% CL. 
Their approximate analysis neglects baryons (this has also been noted recently by Beca et al. 2003) and is based on a linear 
treatment of perturbations until the present time. However, the nonlinear clustering of GCG could be rather nontrivial because 
of the time dependence of its Jeans length and deserves more study. 

In Fig. 8 we compare the predicted power spectra - total (solid line), GCG (dashed line) and baryon (dotted line) - for 
UDM model with the APM (Padilla & Baugh 2003) and 2dF (Colless et al. 2001) data when a varies. The power spectra 
are obtained through the transfer functions of the modified CMBFAST code (CF), i.e. by using the linear evolution. We do 
not show as target the CDM power spectrum because almost indistinguishable from the baryon one. The power spectra are 
normalized at k — 0.01 h/Mpc as in Sandvik et al. (2002) and the results in Fig. 8 can be directly compared with those of 
Fig. 1 in Sandvik et al.. It is very interesting to note the role of baryons, which were neglected by Sandvik et al.: baryons 
keep on clustering on all scales, at all time after decoupling. Even if baryons are a subdominant component (Qb = 0.05), 
they are very important when one considers the total matter power spectrum Pr(k) = f^Pfc(fc) + Q x Px(k). The inclusion of 
baryons smooths out the oscillations of the GCG component in the total matter power spectrum, in particular for a > 10 -3 : 
when comparing to the LSS data, only the shape is relevant and not the features of the GCG component. We also note that 
for a > 10~ 2 the appearance of a finite Jeans length and/or of a suppression of DE component maybe very interesting in 
connection with the bump present in APM data. Therefore, it appears that an intermediate value of a (0.01 < a < 0.2) may 
fit both CMB and LSS data. 

On concluding, in UDM models the comparison of the total matter power spectrum with LSS changes by including 
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Figure 6. SNIa likelihood in the plane wx ,ce, with a Gaussian prior Q m = 0.3 ± 0.1. Confidence regions at 65%, 95% and 99% from 
dark to light gray. 

baryons, even within the linear approximation. In the context of DE models, this effect is even more important since CDM 
perturbations clustering on all scales are also present. 



8 CONCLUSIONS 

We have analyzed quantitatively the observational effects of a GCG as a DE candidate. We have restricted our analysis to a 
flat universe and to purely adiabatic pressure perturbations. 

So far, most authors studied only the comparison of the Hubble law predicted in presence of a GCG with the SNIa 
(Fabris et al. 2002a; Avelino et al. 2002; Makler et al. 2003). However, as also our Figs. 6-7 show, SNIa data constrain GCG 
models very weakly. For instance, the Hubble law is almost insensitive to a, as we see in Fig. 6. Instead, a is the maximum 
value which the sound speed can reach and it is related to the time at which the Jeans instability disappears and therefore 
fluctuations are very sensitive to this quantity: this results in a strong dependence of the CMB spectrum on a, as we see from 
Fig. 2. 

The combined result of CMB and SNIa data is presented in Fig. 7, which leads to a < 0.2 and wx < —0.8 at 95% c.l. 
The bound imposed by WMAP data on a is 5-10 stronger than the previous data (Bean & Dore 2003), and the Chaplygin 
gas (q = 1) is ruled out as a DE candidate at more than the 99.99% c.l. 

The possibility of a unified dark model with a GCG (corresponding to fi c = 0) seems disfavoured by the latest CMB 
anisotropy data in comparison to the GCG playing the role of a DE model, although the UDM on its own is not a bad fit to 
the data. The x 2 statistics is 987/893 for a UDM model, whereas is 979/892 for a DE model. Our result is in contrast with 
the analysis of the peaks position recently performed by Bento, Bertolami and Sen (2003) , in which a different region in the 
(a,wx) plane is allowed. This disagreement may come from an inaccurate analytic estimate of the peaks position (already 
pointed out by CF). 

We have also addressed the problem of the comparison with the LSS data. We have compared results on the matter power 
spectrum with 2dF and APM data, within the linear approximation. By considering the worst case, i.e. a UDM model, we 
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Figure 7. Adding the SNIa constraint. The thick gray curves are the confidence regions for the SNIa; the dotted curves the confidence 
regions for the CMB; the gray regions represent the combined likelihood (65% and 95% c.l.). As it can be seen, the SNIa constraints add 
little information. 

find that LSS data may be complementary to CMB data. The main question about the validity of linear treatment on such 
short scales however remains. 

On concluding, standard quintessence models seem the most economic variation of DE models with respect to ACDM. 
Differentiating among DE models by studying the sound speed of its perturbations seems a promising avenue and this has 
already been shown in two different class of models (CF; DeDeo et al. 2003). We have shown that GCG models are models 
with high predictive power when compared with CMB anisotropics. Current and future * high precision CMB anisotropy data 
can finally constrain and/or rule out physical scenarios based on different models of DE. 
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